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Abstract 



We study the A0f +1 kink soliton and the zero-mode contribution to the 
kink soliton mass in regions beyond the semiclassical regime. The calcula- 
tions are done in the non-trivial scaling region and where appropriate the 
results are compared with the continuum, semiclassical values. We show, as 
a function of parameter space, where the zero-mode contributions become 
significant. 
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Solitons are non-dispersive localised packets of energy moving uniformly, and resembling 
extended particles. The elementary particles in nature are also localised packets of energy, 
being described by some quantum field theory. Because of these features, the soliton might 
appear as the ideal mathematical structure for the description of a particle. When it was re- 
alised that many nonlinear field theories used to describe elementary particles, also had soliton 
solutions and that these solutions might correspond to particle type excitations, the develop- 
ment of methods for soliton quantisation became important. The quantisation of solitons is 
usually done by performing an expansion in powers of H (loop expansion) such that the classical 
soliton solution appears as the term of leading order in the expansion and terms of higher order 
represent the quantum effects. In the mid 1970's there appeared a number of works which 
developed the semiclassical expansion in the quantum field theory. In this period, there were 
schemes being constructed to quantise these solitons. The correspondence between classical 
soliton solutions and extended-particle states of the quantised theory were established |],[|||. 
Collective coordinate methods were used to deal with the so-called "zero-mode problem" f6)-§. 
This problem is a manifestation of the translational invariance of the theory, broken by the in- 
troduction of the soliton. Field oscillations around this classical solution contain zero frequency 
modes, describing displacements of the soliton. 

Here we study the mass of the simplest topological soliton, that is the \4>\+i kink, using 
lattice Monte Carlo techniques. The dynamics of this model are governed by a Euclidean 
Lagrangian density 

£ = \{d^) 2 -^ 2 + \d>\ (1) 

where \l is a bare parameter and A is the bare coupling constant. For a free scalar field theory 
A — > and — [i 2 — > m 2 , where m would then be the bare mass of the <fi. In the classical theory 
fi 2 > corresponds to the onset of spontaneous symmetry breaking. There are two types of 
static solutions, both being static solutions, to the equation of motion; the trivial solutions 



and the topological solutions 



0k,ak = ±-^tanh[-^=x], (3) 

where k and ak label the topological solutions with the kink and antikink correspond to the 
positive and negative sign, respectively. The semiclassical regime corresponds to large /. 

The kink provides an example of a particle with an internal structure, distributed over a 
finite volume rather than concentrated at one point. It possesses a nonzero, conserved quantity 
called the topological charge Q which is defined as 

Q = W*)L-*(*)LJ 
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and which leads to stability of the kink solution. The classical kink mass M c i, is defined to be 
the energy of the static soliton and is given by 



The vacuum and kink solutions can be quantised by path integral quantisation or by construc- 
tion of a tower of approximate harmonic oscillator states around the classical solution 0, where 
either = <p , <p k or a k. In a finite box with the length L the soliton mass with a one loop 
quantum correction becomes [Hj: 



M sol = E kiDk - E vac = M cl + J2 \"l ~ E \& + ( 5 ) 
where uj n and are the eigenvalues of the following equation 

d 2 



dx 2 



+ (3A0 2 - 2/i 2 )_ 



r]i(x) = e^i, (6) 



with = O and = 0^ or a k for u n and £ ra respectively. An important point to note is that, 
due to translational invariance, one of the cj's is zero, i.e, ujq = 0. 

As L is set to infinity then the sums are replaced by integrals and one ends up with a 
logarithmically divergent integral and renormalisation is required to render the kink mass finite. 
Then one arrives at the mass of the continuum kink with one loop corrections [[|: 

M^^ + M^f-^W). (7) 

The zero eigenvalue to = is referred to as the zero-mode and has well known physical conse- 
quences. These modes always exist when one quantises a theory with a translationally invariant 
Lagrangian about a solution that is not translational invariant. The physical consequence of the 
existence of the zero-mode is the center of mass motion of the kink. In Eq. (0) the zero-mode 
contribution to the kink mass is omitted. That effectively means that the mass of the quantum 
kink particle is assumed to be the same as the kink energy. The semiclassical treatment of these 
zero-modes is done in variety of ways however, we will not discuss these further here, 

since it is not necessary for our purposes. It is also important to mention that the semiclassical 
quantisation of the discrete lattice version of the lagrangian density given by Eq. (|1|) in a finite 
size box is also complicated by the zero-mode problem. 

Kink on the lattice 

The A0 4 action on a 2-d lattice can be written as: 

s = -J2 <^<W + E( 2 - f )ti + -A (8) 
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where we have defined the dimensionless quantities p = fia and A = Aa 2 with a being the 
lattice spacing. In addition n = (n 1; n 2 ) is a 2-dimensional vector labeling the lattice sites and 
fi is a unit vector in the temporal or spatial direction (not to be confused with our parameter 
p 2 = — m 2 ). We have also denoted the field on the neighboring site of n in the direction of \x 

by 0n, M - 

This model exhibits two phases. In some regions of the phase space ((f)) = and these are 
the symmetric (unbroken symmetry) regions, whereas in other regions spontaneous symmetry 
breaking occurs and ((f)) ^ 0. Classically, for positive values of m 2 = —p 2 one always has 
((f)) = and for negative m 2 (i.e., positive p 2 ), where spontaneous symmetry breaking occurs, 
(4>) 0. In this regime the second order critical line which separates the two phases is the line 
corresponding to m 2 = 0. Beyond the classical limit the phase space structure changes. There 
are still two phases and there is still a second order phase critical line separating these phases, 
however, the position of the critical line changes and due to quantum fluctuations washing out 
"shallow" spontaneous symmetry breaking it occurs at a finite negative fh 2 in general. 

In order to determine the critical line, we choose several values of m 2 < located in the 
broken symmetry sector where (<fi) ^ 0. For each value of m 2 , A can be increased until ((f)) = 
and, thus, the critical parameters, (m 2 ,A c ), can be found. Of course there is no second order 
phase transition on a finite lattice, however, by a second order phase transition here we mean 
that the correlation length is much larger than the lattice dimensions. The critical line is shown 
in Fig. 1. We have also shown in this figure the one loop prediction for the critical line using 
the light-front formulation ||15|| . 

The authors of Ref. [16| proposed two methods of calculating the kink mass on the lattice. 
Here we use one of these methods which uses a local parameter and is much less susceptible 
to finite size effects. In this method one imposes an anti-periodic spatial boundary condition, 
giving rise to a topological excitation with non-zero topological charge. Since the kink has the 
lowest energy in the topological sector this topological excitation corresponds to the kink. It is 



shown that for a a fixed m, one has fl6|j : 

M 8ol 0) = f 0tl(fi?) = \ \1 dff\ \(S a ) - (S p )} , (9) 

where M so \ is the soliton mass, /3 — 1/A, /3 C is the inverse of the dimensionless critical bare 
coupling A c , T is the length of the lattice in the temporal direction, (S p ) and (S a ) are the mean 
action of the system with a periodic and anti-periodic spatial boundary condition, respectively. 

As we mentioned earlier in the semiclassical quantisation one encounters the zero-mode 
problem with its physical consequences being the center of mass motion of the kink. An 
interesting question is whether this problem persists beyond the semiclassical regime. To answer 
whether the zero-mode problem persists beyond the semiclassical regime, one can examine one 
of the consequences of the existence of a zero-mode, that is, the kink displacements. On a 
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lattice with anti-periodic boundary condition in the special direction, we set m 2 = — fi 2 = — 1 
and A = 4 giving / = 0.5, which corresponds to a region beyond the semiclassical regime. Then 
for an arbitrary time slice we calculated (<f) n ) for each site for a number of configurations and 
an average over these configurations was calculated. We have shown the results in Fig. 2. As 
this figure suggests, the movement of the kink due to the translational mode still persists even 
beyond the semiclassical regime. We repeated the same procedure for different time slices and 
different number of sampled configurations and these results showed the same behaviour. 

In order to treat the problem, we imposed an additional constraint, that is <fr(M) = with 
M = (xq,N/2) for all times Xq, fixing the center of the kink to the center of the lattice. For 
the same time slices as in the previous case and for a number of constrained configurations, 
the value of (0 n ) for each site along with the classical configuration are included in Fig. 2. The 
constrained configuration resembles the kink solution with its center located at the center of 
the lattice. 

Both in numerical MC studies and the analytical calculations it is important to find the 
renormalisation group trajectories (RGT). Along these curves and close to an infrared (IR) 
fixed point the physics described by the lattice regularised quantum field is constant and only 
the value of the cut-off (lattice spacing) is changing. This region is called the scaling region. 
The A0 4 theory in two dimension is an interacting theory. That is, in addition to having a 
Gaussian fixed point at m 2 = A = where the renormalised coupling A r vanishes, it has other 
infrared fixed points at which A r is non-zero. The best candidates for these fixed points are 
along the critical line where a second order phase transition occurs. In this model the only 
non trivial critical region is along the transition line shown in Fig. 1. In the vicinity of the 



fixed point the vertex functions strongly scale flT] and one expects that close to the critical 
line, there should be segments of phase space where the ratios of dimensionless vertex functions 
remain almost constant, giving the scaling region. 

In our calculations it is important to find the scaling region corresponding to a non-trivial 
IR fixed point. That is, one should try to find trajectories away from the Gaussian fixed 
point. On trivial fixed point trajectories, even though spontaneous symmetry breaking can 
still occur and hence a kink solution can exist, the vacuum is governed by a free field. We 
used R(m r ,\ r ) = m 2 /X r as a dimensionless parameter for probing the scaling region. This 



parameter R, can be calculated accurately using the effective potential method [18]. The 



scaling region corresponds to regions where R is approximately constant. Our entire calculations 
are performed within the scaling region so that they can be legitimately compared with the 
continuum-semiclassical predictions for the soliton mass, given by Eq. (|?|). 

For a fixed value of m 2 = — 1 we found a range of values of A in the broken symmetry 
sector, that is 0.2 < A < 0.8, for which the values of R(m r ,X r ) were approximately constant, 
determining a segment of the scaling region corresponding to m 2 = — 1. Then for some values 
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of A within this region we calculated AS = (S a ) — (S p ). In calculations of (S a ) we constrained 
the center of the kink to the center of the lattice. We have plotted AS/T versus f3 along with 
its classical value in Fig. 3. Having known AS/T , Eq. (|]) can be used to calculate the soliton 
mass. We also calculated the soliton mass with and without imposing a constraint on the 
center of the kink. The comparison of these results with each other and with the classical and 
semiclassical values is shown in Fig. 4. As is evident, in this case the zero-mode contribution 
is within statistical uncertainties. The statistical uncertainties, as one expects, increase as one 
approaches the critical line which makes the detection of the zero-mode contribution to soliton 
mass difficult. Also, it is interesting to note that the Monte Carlo results for the soliton mass 
are less than the classical mass but larger than one loop semiclassical predictions. 

Next, for m 2 = —2.2, we repeated the same procedure and calculated the soliton mass for 
a number of bare couplings within the scaling region. The calculations were done with and 
without a constraint on the center of the kink and the results and their comparison with the 
classical and semiclassical predictions are shown in Fig. 5. Our results are consistent with 
Ref. [Kj. There are two important observations. First, unlike the previous case the MC 



calculated soliton masses are lower than the semiclassical values which is due to the higher 
order corrections. The other important point is that, close to the critical line it appears to be 
clear that there is a positive contribution to the soliton mass due to the zero-mode. 

Finally, for m 2 = —4 we calculated the soliton mass, with and without a constraint on the 
configurations. The results are shown in Fig. 6 and again the zero-mode contribution to the 
soliton mass appears to be positive. 

All our calculations were done on a 48 x 48 lattice. As one approaches the critical line the 
correlation length increases and the finite size effects become more significant. However, since 
the Monte Carlo calculation of masses are based on a local parameter fl(ft) the finite size effects 
are smaller than one might, in general, expect fl9 |. In our calculations we kept the correlation 
length below a half of the lattice length. 

In conclusion we would like to mention that in addition to a straightforward elimination of 
the zero-mode, the imposition of a constraint on the center of the kink resulted also in more 
stable configurations and a reduction of the statistical uncertainties on (S a ) and consequently on 
M so i. These instabilities are more significant close to the critical line where the field fluctuations 
are larger. In general we found that the statistical uncertainties on (S a ) where much larger than 
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FIGURES 
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FIG. 1. The plot of the transition line between the broken sector and unbroken sector using Monte 
Carlo methods (diamond) and the light-front perturbative predictions (dashed line). We used the 
symbol r = m 2 in the above figure. 
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FIG. 2. The comparative plot of averges of fields on a time slice with / = 0.5 and m 2 = — 1 for 
a unconstrained lattice (dashed line) and constrained lattice (dotted line). The constraint fixes the 
center of the kink to the center of the lattice for all time slices. The solid line is the classical solution. 



8 



1 .0 



0.! 



< 



0.6 



0.4 



0.2 



0.0 







2 



$ $ * 







FIG. 3. The plot of AS/T versus = l/\ with m 2 
classical value for AS/T 
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FIG. 4. Plot of soliton mass versus (3 with fh 2 = —1. The Monte Carlo results are compared with 
the classical (solid line) and semiclassical results (dashed line). 
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FIG. 5. Plot of soliton mass versus (3 with to 2 = —2.2. The Monte Carlo results are compared 
with the classical (solid line) and semiclassical results (dashed line). 
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FIG. 6. Plot of soliton mass versus (5 with m 2 = —4. The Monte Carlo results are compared with 
the classical (solid line) and semiclassical results (dashed line). 
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